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Testing the spin-cutoff parameterization with shell-model calculations 

William M. Spinella and Calvin W. Johnson 
Department of Physics, San Diego State University, 
5500 Campanile Drive, San Diego, CA 92182-1233 

The nuclear level density, an important input to Hauser-Feshbach calculations, depends not only 
on excitation energy but also on angular momentum J. The J-dependence of the level density 
at fixed excitation energy E x is usually parameterized via the spin-cutoff factor a. We carefully 
test the statistical accuracy of this parameterization for a large number of spectra computed using 
semi-realistic interactions in the interacting shell model, with a nonlinear least-squares fit of a and 
finding the error bar in a. The spin-cutoff parameterization works well as long as there are enough 
states to be statistical. In turn, the spin-cutoff factor can be related to the average value of J 2 at 
a fixed excitation energy, and we briefly investigate extracting {J 2 {E X )) from a thermal calculation 
such as one might do via Monte Carlo. 

PACS numbers: 
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An important input into compound nuclear reactions such as statistical neutron capture is the level density [I]. 
The level density is deceptively simple: the number of levels per unit of excitation energy. But when thousands or 
millions of states or more are involved, as in the statistical regime, experimental measurement becomes challenging 
and theoretical calculations are also difficult if not intractable. Because of this the level density is one of the most 
uncertain inputs into reaction calculations [2j. 

(Note: one should be careful to distinguish between the level density, which does not include the 2 J 4- 1 degeneracy 
^ . in M , and the state density which does include the 2 J + 1 degeneracy in M.) 

For calculation of statistical neutron capture rates one needs not only the density of levels p as a function of excitation 
energy E x but also as a function of parity and angular momentum J. While there are different parameterizations for 
the level density p{E x ), such as the back-shifted Fermi-gas [3|-[5[ or constant temperature [4J], for the dependence on 
^. \ J the parameterization [fj 

\o'- 

j5 : p(E x ,J) = p(E x )f(J), /(J) = ^±i cxp (-J(J + l)/2a 2 ), (1) 

I* ■ is universally used, where a 2 is called the spin-cutoff factor. (This normalization is for the level density; for the state 
(~ } density the spin-cutoff parameterization has a different normalization.) 

A number of other works have investigated the spin distribution of levels, both experimentally [3] and theoretically 
t-H ' [8|, |9| . While these papers have compared the actual spin distribution against the spin-cutoff parameterization, in this 
^ . paper we focus on how well does the spin- cutoff parameterization actually describe the spin distribution, in particular 
by computing the error bars on the spin-cutoff factor. 

Towards this end, we have generated a number of theoretical nuclear spectra, as described in more detail in Section 
HT1 We then put the data into energy bins and, in each energy bin, did a nonlinear least-squares fit for a using Eq. ([1} , 
including deriving the error bars on a. The fitting methodology and results we give in Section [IIII 

In general we find the unsuprising but gratifying result that the spin-cutoff parameterization does a good job of 
describing the spin distribution, and that the error bars on a decrease as the number of levels in a bin increase. If 
Eq. ([1} were wrong, the error bars would not decrease systematically. 

Because Eq. ([1]) has only the one parameter a, it is tempting to try to extract the spin-cutoff factor directly rather 
than fitting it, for example relating a 2 to the average value of J(J + 1) at a given excitation energy. This would be 
useful in particular for Monte Carlo calculations such as the so-called shell- model Monte Carlo (SMMC); although a 
J-projected version of SMMC has been developed, it is very computationally intensive. By contrast, it is relatively 
easy to compute thermally-averaged expectation values of operators such as J 2 in SMMC. Because one needs the 
average as a function of excitation energy and not temperature, one needs to extract the energy-averaged (J(J+1)(E)) 
from thermally-weighted values, in the same way one inverts the Laplace transform to extract the state density from 
the thermal partition function. This we discuss in Section [IVI 
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II. SHELL-MODEL FRAMEWORK 



We work in the framework of the configuration- interaction shell model [1CH12]]. Here one defines a finite single- 
particle space and has as input single-particle energies and two-body matrix elements. We use semi-realistic/semi- 
phcnomenological matrix elements, meaning they start from realistic fits to scattering data, are then renormalized, 
and finally tuned to reproduce binding energies and spectra (see [lCj, LL3| for details of methodology) . 

All model spaces assume some inert core and valence particles. The three model spaces we work in and their 
interactions are: 

• sd, or the lsT/rOdz /2-0e? 5 / 2 valence space, assuming an inert 16 core; the interaction is the universal sd-interaction 
'B,' or USDB |i|; 

• pf, or the I-P1/2-I.P3/2-O/5/2-O./V/2 valence space, assuming an inert 40 Ca core; the interaction is the monopole- 
modificd Kuo-Brown G-matrix interaction version 3G, or KB3G |14l.[l5|: 

• and finally the ps-d^/2, or 0p3/2-0j>i/2-0<i5/2- lsi/2 valence space, assuming an inert 4 He core; the interation is a 
hybrid of Cohen-Kurath (CK) matrix elements in the Op shell[la|, the older universal sd interaction of Wildcnthal [13] 
in the Od 5 / 2 -ls 1 / 2 space, and the Millener-Kurath (MK) p-sd cross-shell matrix elements [18j|. We leave out the 0d 3 / 2 
orbit to make calculations tractable. Within the p and sd spaces we use the original spacing of the single-particle 
energies for the CK and Wildenthal interactions, respectively, but then shift the sd single-particle energies relative to 
the p-shell single particle energies to place the first 3~ state at approximately 6.1 MeV above the ground state. The 
rest of the spectrum, in particular the first excited + state, is not very good, but the idea is to have a non-trivial 
model, not exact reproduction of the spectrum. This model space and interaction allows us to consider model nuclei 
with both parities and to investigate parity-mixing in the HF state. 

Finally, we computed the many-body spectra using the BIGSTICK configuration-interaction diagonalization code 
[13 ] . Table U lists the nuclides and the single-particle spaces we computed. 



III. RESULTS 



To analyze the spin distribution of shell-model spectra, we binned the data into 1 MeV bins (our results were 
insensitive to the size of the bins, as we discuss at the end of this section), and did a non-linear least-squares fit to the 
Ericson formula using the MP FIT code [2(J, [2l[ , an implementation of the Levenberg-Marquardt algorithm [22| . For 
the purpose of the fit we estimate the error of each J-projected density pj in an energy bin to be the standard counting 
statistical error, Apj = v / pJ (bins with pj = are assigned an error of 1). The error bar on the least-squares-fit 
spin-cutoff parameter is just the standard square root of the diagonal of the covariance matrix, which in our case was 
a trivial lxl matrix. To reduce overfitting, when our \ 2 per degree of freedom, or reduced \ 2 , 
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fell below 1, we scaled the errors until Xrcd = 1- 

We analyzed spectra from a wide range of nuclides in three different model spaces, and considered even-even, 
odd-odd, and odd- A nuclei. We will show representative plots from a number of different calculations. 

We begin in the sd-shell, where we consider 26 A1, 26 Mg, and 33 P in Figs. Q] [21 and El respectively, showing the 
distribution of states as a function of J. Because we extract error bars on our spin-cutoff factors, we display not only 
the least-squares fit but also the one-standard-deviation envelope (the dark gray band in the figures). The points 
represent our binned CI shell- model data and the error bars represent the statistical error y/pj on the counts in a 1.0 
MeV bin. Not only do the error bars on the CI data get smaller, as the number of counts in the bins increase, our 
error on a decreases as well. 
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FIG. 1: Al in the sd shell. The error bars are standard JJTj statistical error in counting the number of states in a 



given bin of energy 1.0 MeV and fixed J. 



In Fig. 2]we show the evolution of the spin-cutoff factor as a function of for these nuclides as well as several additional 
sd-shell nuclides. Fig. [5] does the same for pf shell nuclides, while Fig. [6] does for both positive and negative parities 
in the p-sd 5 /2 space. The error bars are larger for the last case because the statistics is smaller, and the values of 
the level density in all cases are no doubt an underestimate because intruders are left out. We remind the reader 
our goal is not to determine the 'true' value of the spin-cutoff factor but to compare the statistical goodness against 
semi-realistic calculations. 

FigureslHEJ and[5]show that as excitation energy increases the distribution in angular momentum is better described 
by the spin cut-off parameterization Eq. (fT]), consequently resulting in the reduction of the error on the extracted 
spin-cutoff factor. The pattern of improvement in the fit is observed in all analyzed nuclides, and can be seen by 
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FIG. 2: 2b Mg in the sd shell 



examining the root mean square error (RMSE) between the CI data and the fit: 
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(3) 



The degree of improvement in (pj//o)rmse appears to be independent of any structural differences in the sd and 
pf shell nuclides, as illustrated in Fig. [7] The cross-shell nuclides in the p-sd 5 / 2 space, Fig. [5] exhibit a similar 
improvement in (pj/p)bmse, though because the levels are split between parities the level density is lower than that 
of the sd and pf shell nuclides, and thus the improvement occurs at higher excitation energies than with the latter. 

We also investigated the dependence upon the bin size. While the above results were for standard bin sizes of 1.0 
MeV, we also used bins of 0.5 and 2.0 MeV. The best-fit values were insensitive to the bin size, as illustrated in Fig. |H1 




0.30 



0.25 



0.30r 




0.30 




2.8 ± 0.1 



15.0 < E x (MeV) < 16.0 



0.30 



0.25 



0.20 



3.17 ± 0.07 



16.0 < E x (MeV) < 17.0 




FIG. 3: 33 P in the sd shell 



although the error bars depended upon the bin size, which is primarily an effect of the number of states in a bin. 

We therefore conclude that the spin-cutoff factor is a statistically effective parameterization of the J-dependcncc 
of the level density. 



IV. MONTE CARLO CALCULATIONS 



We attained our level densities above through laborious diagonalization of shell-model Hamiltonians, but such 
calculations are not practical for large model spaces. Therefore we consider other techniques to compute the level 
density and the spin-cutoff factor. 

Most theoretical calculations of the density of states start from thermodynamics, as the partition function is the 
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FIG. 4: Fit spin-cutoff factors for selected sd-shell nuclides 



Laplace transform of the density of states: 



Z(P) 



tve-P" = f e- 0E p(E)dE = J2 e M-/3E l 



(4) 



where the trace and the sum are over many-body states and one can interpret /3 as inverse temperature. The density 
of states itself can be written as p{E) — '^ ji 5{E — Ei). The Laplace transform can be inverted with reasonable 
reliability through the saddle-point approximation [5j, and so now one has the problem of computing the nuclear 
partition function. 

One avenue to compute the many-body partition function is to use auxiliary-field path integrals to evaluate 
exp(— j3H). While such an approach has been successful, it requires a nontrivial amount of computing time; further- 
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FIG. 5: Fit spin-cutoff factors for selected p/-shell nuclides 



more because one must evaluate the resulting integrals through Monte Carlo sampling, one has to use a well-posed 
interaction that is free or mostly free of the sign problem [24| . 

The spin-cutoff factor can be related directly to the average value of J(J + 1) in an energy bin. One can either 
assume the nucleon spins add together randomly [25| or simply replace the discrete Ericson function ([l) by a continuous 
function of J and integrating. Then one finds 



2a 2 « ( J( J + 1)> + - 



(5) 



for the level density, or 3cr w ( J( J + 1)) if one is using the state density as in [25[. Using our level densities, we found 




. 


"1 


" 


" 


1 


1 




"■■ 


12 c - 





10 15 20 25 30 
Excitation Energy (MeV) 



10 15 20 

Excitation Energy (MeV) 



25 



.I'l'l i 



1(1 fj'sf*'*" 1 



14 



N" 



10 15 20 

Excitation Energy (MeV) 



25 30 



4.0 
3.5 
3.0 
2.5 
b2.0 
1.5 
1.0 
0.5 
0.0, 



' 




\ 


J I 1 * ' ' - 


' 1} 


I ' 


17 





10 15 20 

Excitation Energy (MeV) 



25 



4. Op 
3.5 
3.0 
2.5 
b 2.0 
1.5 
1.0 
0.5* 
O.oi- 



IK'/"' 



\ 



20 



Ne~ 



10 15 20 

Excitation Energy (MeV) 



25 30 



\ I 



i i 



i *i 



i i 



21) 



Ne _ 



10 15 20 

Excitation Energy (MeV) 



25 



FIG. 6: Fit spin-cutoff factors for selected p-sd,5/2-she\l nuclides, for both positive parity (left-hand figures) and 

negative parity (right-hand figures) 



this approximation worked well, as illustrated in Fig. 1101 

The authors of [251 ] relate the spin-cutoff factor to the moment of inertia which then they compute in rather 
complicated fashion. However, it is straightforward to compute the thermal average of the expectation value, 
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(J 2 W) =' L g , m i=Z(/3)- 1 ^exp(-/3^)J,:(J t + l). 
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FIG. 7: (p,//p)rmse for all analyzed sd and pf shell nuclides organized by structure. 



This is just the Laplace transform of the sum of J( J + 1) in an energy bin, that is if we define the average as 



(J(J + 1)(E)) = 



P(E) 



(7) 



(we deliberately use the notation J 2 for the thermal average and J(J + 1) for the energy average, to help distinguish 
the two notationally) then 

Z(/3)(J 2 (/3)) = [ e-' iE p{E){J(J + l){E))dE. 



(8) 

In the same way one can invert the Laplace transform via the method of steepest descents, one can also extract 
(J (J + 1)(E)), or rather, extract p(E){J(J + l)(E)) and then divide by the level density, as follows. First, find the 
saddle point, which yields the effective temperature for a given target excitation energy: 
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For a given effective temperature T c s = l//3 e ff, the energy E is different from that for the saddle point condition for 
the level density, because of the additional term. 

One must also compute the effective specific heat, which is the second derivative: 
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The final result of the method of steepest descents is 



p(E)(J(J+l)(E)) = 



Z(/3 cff )(J 2 (/3 cff ))cxp(/3 cfF £) 



(10) 



(11) 



The level density itself is computed using the exact same kind of inversion, only without the (J 2 (/3)) term in the 
derivatives. Although we did not carry out any SMMC calculations, we modeled such a calculation by computing 
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FIG. 8: (pj/p)rmse for all analyzed psds cross-shell nuclides. 



an exact and complete spectrum, in this case for 44 Ti in the pf shell using the KB3G interaction, which has an M- 
scheme dimension of 4000; from that spectrum we constructed the partition function as well as the thermal average 
of J(J + 1). We then inverted using Eq. (fTTj) and got a good value for a 2 , as show in Fig. [TTJ 

Although a practical J-projection scheme for level densities has been implemented [25J, by computing an M- 
projected density and taking pj — pu=j~ pm=j-i, simply computing ( J 2 (/3) is significantly computationally cheaper, 
and appears to be effective in arriving at the spin-cutoff factor. 

V. CONCLUSIONS AND ACKNOWLEDGEMENTS 

We have carried out a detailed investigation into the spin-cutoff parameterization Eq. (fTJ) of the J-dependence of 
the level density, with particular attention paid to the error bars in the fit to a. We confirmed that the spin-cutoff 
parameterization is a good one, and that as the statistics improve the error bars on a decrease. There was no 
qualitative difference in behavior for different nuclides. 

We also showed that a simple energy average, (J {J + 1)(E X )) gives a good value for the spin-cutoff parameter 
and that one can extract the energy average from a thermal average (J 2 (/3)) which one might get naturally out 
of a quantum Monte Carlo calculation. Such an approach would be computationally more efficient than directly 
calculating J-projected densities, although the latter are achieveable. 

The U.S. Department of Energy supported this investigation through grant DE-FG02-96ER40985. 
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